Eigenvalue Spectrum of a Dirac Particle in Static and 
Spherical Complex Potential 



Khaled Saaidi^ 

Department of Science, University of Kurdistan, Pasdaran Ave., Sanandaj, Iran 

Institute for Studies in Theoretical Physics and Mathematics, P.O.Box, 19395-5531, Tehran, Iran 



Keywords: Complex potential, non-Hermitian, Dirac equation, energy spectrum. 



Abstract 

It has been observed that a quantum theory need not to be Hermitian to have a 
real spectrum. We study the non-Hermitian relativistic quantum theories for many com- 
plex potentials, and we obtain the real relativistic energy eigenvalues and corresponding 
eigenfunctions of a Dirac charged particle in complex static and spherically symmetric po- 
tentials. Complex Dirac-Eckart, complex Dirac-Rosen-Morse II, complex Dirac-Scarf and 
complex Dirac-Poschl- Teller potential are investigated. 
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1 Introduction 



The first interest study in the non-Hermitian quantum theory date back to an old paper by 
Caliceti et al [1] . The imaginary cubic oscillator problem in the context of perturbation theory 
has been studied in [1]. The energy spectrum of that model is real and discrete. It shows that 
one may construct many new Hamiltonians that have real spectrum, although, their Hamilto- 
nians are not Hermitian. The key idea of the new formalism (non-Hermitian quantum theory) 
lies in the empirical observation that the existence of the real spectrum need not to necessarily 
be attributed to the Hermiticity of the Hamiltonian. Such non-Hermitian formalism, for the 
context of Schrodinger Hamiltonian have been studied for several models in [1-25] . The analysis 
of the related purely real spectra of energies has been performed with different techniques. For 
example, resummations of divergent perturbation series[l], delta expansions [2], WKB method 
[6], functional analysis [8], pseudo-Hermitian methods [17, 21, 23], and complex Lee algebra 
[18-20]. Also, some explicit studies of the Hermitian and non-Hermitian Hamiltonians have 
performed in the context of Dirac Hamiltonian. For example, the solution of ordinary (Her- 
mitian) Dirac equation for Coulomb potential including its relativistic bound state spectrum 
and wave fimction was investigated in [22,23]. Also by adding off-diagonal real linear radial 
term to the ordinary Dirac operator, the relativistic Dirac equation with oscillator potential 
has been introduced [24, 25] then the energy spectrum of corresponding eigenfunctions have 
been obtained. The ordinary (Hermitian) Dirac equation for a charged particle in static elec- 
tromagnetic field, is studied for Morse potential [29]. The non-Hermitian formalism, for the 
context of Dirac Hamiltonian, have been studied for three dimensional complex Dirac-Morse 
and complex Dirac-Coulomb potentials in [30, 31]. 

In this paper, we consider the non-Hermitian Dirac Hamiltonian for several complex po- 
tentials. We study a charged particle in static and spherically symmetric four vector complex 
potentials. By applying a unitary transformation to Dirac equation, we obtain the second order 
Schrodinger like equation, therefore comparison with well-known non-relativistic problems is 
transparent. Using correspondence between parameters of the two problems (the Schrodinger 
equation and the Schrodinger like equation which is obtained after applying the unitary trans- 
formation to Dirac equation for a potential ) we obtain the bound states spectrum and corre- 
sponding eigenfunctions. 

The scheme of this article is as following. In sec. 2, we study the non-Hermitian version of 
Dirac equation for a charged particle with static and spherically symmetric potential, and by 
applying a unitary transformation we obtain the proper gauge fixing condition and Schrodinger 
like differential equation. In sec. 3, we discus the Dirac equation for complex Dirac-Eckart 
potential and obtain the real energy spectrum and corresponding eigenfunctions. In sec. 4, we 
consider the Dirac equation for a complex Dirac-Rosen-Morse II potential then we obtain the 
real eigenvalue and its wave function. In sec. 5, and sec. 6, we studies the complex Dirac-Scarf 
and Dirac-Poschl- Teller potentials respectively. We obtain the relativistic energy spectrum and 
correspondening wave function for the upper component of spinors. 

2 Preliminaries 

The Hamiltonian of a Dirac particle for a complex electromagnetic field is (c = ?i = 1) 

H = a.{p - eA) + i3m + eV, (1) 
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where the Dirac matrices d, /3 have their usual meaning, and setting Ao equal to V. In (1), A 
and V arc the vector and scalar complex field respectively, where. A* ^ A and V* ^ V . Then 
the Dirac Hamiltonian (1) is not Hermitian. It is well known that the local gauge symmetry in 
quantum electrodynamic implies an invariance under the transformation as: 



{y,A)^{V',A') 



(2) 



Here A(t, "f*') is a complex scalar field. Suppose that the charge distribution is static with 

spherical symmetry, so the gauge invariance implies that V' = V and A' = fA{r), where r is 
the radial unit vector[29]. One can denoted the correspondence wave function of (1) as: 



* = 



(3) 



In this case one can obtain 



{m + eV - Er)^ = 



{eV -m- Er)x = 



a.^ — e{a.f)A{r) 
a."^ + e{a.f)A{r) 



(4) 



Here it's are the three Pauli spin matrices, Er is relativistic energy, and we replaced iea .A{—iea .A) 
in first(sccond) equation of (4) instead of ea.A, respectively. Note that, because of the spherical 
symmetry of the complex field, the angular-momentum operator J and the parity operator, P, 
commute with the Hamiltonian and the two spinors $ and x have opposite parity also. So the 
correspondence wave functions are denoted by 



X = /(r)CTrO_«,^(^?,(p). 



It is seen that 



, Ik 
(a.^)(/(rKn_«,^(^?,V')) 



r r 

1 K 

(^r^-K,n{^r H )f{r) 

r r 



where k is the spin orbit coupling operator as: 

k — a.L 

and we have used from 



hi. 



in which 



-{l + l) = -{3 



l = ij + 



for j=l+^ 
for j = 1-1 



(5) 

(6) 
(7) 

(8) 
(9) 
(10) 
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Therefore by defining m — g{r)/r , U2 = f{r)/r, we obtain the following two component radial 
Dirac equation [32] 



(m + eV — Er)ui(r) 



{dr~--eA{r))u2{r), 



k 



{eV - m - Er)u2{r) = -{dr -\ \- eA{r))ui 



(11) 



Note that, A{r) is a gauge field, which has a symmetry such as (2), therefore, it must be 
fixed. It is seen that fixing this gauge degree of freedom by "^.^ = ^ = is not a suitable 
choice. Remark that in this paper, instead of solving Dirac equation we want to solve the 2*^ 
order differential equation, which is obtained by eliminating one component of equation (11). 
However, for the cases which eV ^ 0, the second order differential equation is not Schrodingcr 
like and, therefore, one can obtain the proper gauge fixing by applying the global unitary 
transformation on two components u\ and ui such as: 



U 



a ib 
ib a 



where a,b G 5R, and + 6^ = 1. By applying (12 ) to the upper component, 
component, of spinor and institute them in (11), we have 



(12) 
and lower 



(m - ErC)(f>'' + 



S 



eV - iSEr - dr 



= 0, 



eV - iSEr + dr 



r - (m + ErC)(j)'- = 0, 



where, S= 2ab, C = — and we have used from a gauge fixing condition as: 

eV=-{eA+-). 



(13) 



(14) 



However, we eliminate the cj)'' component in (13), and obtain the radial differential equation 
for (/)" as: 

'^'^"+Veffr + {m^-K)<l>"=0, (15) 



where 



eff 



{eVy + iEriS"- - C'){eV) - i 



0, 

(^2 _ C2) d{eV) 



(16) 



52 ' ' ' ' " S dr 

Furthermore, from (11), it is easily seen that for the cases which eV = the unitary 
transformation is not necessary. So, we can rewrite Ui and U2 as: 



Ui = 
U2 = 



then, by eliminating 0', one can obtain the Schrodinger like differential equation for radial 
upper component, 0", as: 



Veffr 



[m 



E'rW = 0, 



(17) 
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where 

ye// = (eA(r) + ^)2-|-(e^(r) + ^). (18) 

3 The complex Dirac-Eckart potential 

The complex Eckart potential which holds discrete energy spectrum is [22] 

(j.) = ^) _ 2iB coth(x) , (19) 

smh(a;) 

and the corresponding Schrodinger equation is: 

+ ^.^7 - 2iB coth(x) - e] V>(x) = 0. (20) 
dx'' smh(a;) / 

Here A,Bg^. It is obviously seen that, the Eckart potential is singular in origin, but, one 

can simply avoid their singularities by a local deformation of the integration path. In [22] 
by solving the Schrodinger equation (20 ), the real energy spectrum and corresponding wave 
functions were found as: 

En = jj—:^ -{A-nf, n=l,2,...,nmax<A, (21) 

and 

^^{x) = Ar„(coth(a;) - l)''(coth(a;) + iyPi^^'^-'\coth{x)), (22) 
where 7V„ is a normalization constant, pj^'^'^"^ is Jacobi polynomial ^and 

2/x = {A-n), 

2v = -4^- (23) 
A-n ^ 

Now by defining the complex Dirac-Eckart four vector potential as: 

{eV{r), eA{r)r) = (^iC coth(r), (^ coth(r) -'^)f^, (24) 

and by using (16), we can obtain, 

^e// = ^^^^-^7C0th(r)-r7^ (25) 
coth (r) 

where C e 3i, r? = §(S'^ -C^) and 7 = 2Er{C^ -S'^)C. So that the 2*'' order differential equation 
for radial upper component is: 

''^''"^^ ncoth{r)-{E^ + r]^ -m^)](f>l{r)=0. (26) 



dr"^ coth^(r 



^P^'''"' (a;) is Jacobi polynomial in all of this text 
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Comparing (26) with (20) and then using (21) and (22), one can arrive at the relativistic 
real energy eigenvalue as: 



in which 



rf — {rj — rif']{ri — nf' 



[(1 - 25)77 -n][(l + 25)r?-n] 



(27) 



- r?| < V??^ -m2, (28) 

and the correspondence wave function for ^!>JJ(r) is the same eigenfunction in (22) with a new 
set of parameters as: 

2/i = (jj - n), 

2v = -J2L = _!I (29) 
r]-n 2/1 ^ ^ 

4 The complex Dirac- Rosen-Morse II potential 

The Schrodinger equation for complex Rosen-Morse II potential is: 

-^^ + V^''^ix)^{x)-Ei;{x)=0, (30) 

where 

yCiJM(^) = [{bR + ibif + g2 - l/4]cosch2(a;) - 2g(6^ + i67)cosch(a;) coth(a;). (31) 
One can rewrite (30) as: 

+ V^^''{x)M^) = -{q-n- l/2fMx) = 0, (32) 

which shows that 

E„ = -{q-n-l/2f, (33) 

where 

n = 0,l,2,...n,„„, < (j-1/2, (34) 
and also q > 1/2. Hence, we define the Dirac-Rosen-Morse 11 four vector potential as: 

{eV{r), eA{r)f) = (o, [Ccoth(r) - {tjr + ir]i)cosch{r) - -]r^ . (35) 

Here, C,r]n,T]i G 3?. By making use of (18), we find the effective complex Dirac-Rosen-Morse 
H potential, Vfj:f^^{r), as: 



(??fl+%)'+(C+l/2)'-l/4) 



cosch^(r)-2(C-M/2)(r/K-h%)cosch(r)coth(r)-|-C^. 

(36) 
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Using (36 ), we obtain the following second order differential equation for the upper component, 

+ [(77K + ir,,)2 + (C + 1/2)2 -l/4)]cosch2(r) 



- 2(C + l/2)(77fl + %)cosch(r) coth(r) 



m = iE'^^-m'-Cmr). (37) 



Comparing (37) with Schrodinger equation for complex Rosen-Morse II potential, (31), gives 
the following real relativistic energy spectrum for complex Dirac-Rosen-Morse II potential as: 

Em = Vm^ + C^-iC-nf), (38) 

where n = 0, 1,2, ...jUmax- Here, ^ > 0, and reality of energy spectrum emphasis that Umax 
satisfy 

\nmax - CI < Vm-' + e- (39) 

5 The complex Dirac-Scarf potential 

In [22], the complex Scarf potential and the corresponding Schrodinger equation is given in the 
form 

V^^{x) = [{bR + ibif -q^ + l/4]sech^(a;) - 2q{br + i6/)sech(x) tanh(a;), (40) 



+ [{bR + ibiY -q^ + l/4]sech2(a:) - 2q{br + i6/)sech(x) tanh(x) - E 



i){x) = 0. (41) 



It is well known that the Schrodinger equation of this potential is exactly solvable. It is 
easily seen that this complex potential (40) is not invariant under PT-symmetry, but, for 
the cases which bR = 0, is PT-symmetry invariant (where P denotes the parity operator; 
{Ptp){x) = i'i—x) and T the time reversal operator; {Tip){x) = ip*{x) ). For q > 1/2, the 
associated eigenvalues and eigenfunctions are 

En = -{q-n-l/2f, n = 0,1,2, ...Umax < q - 1/2, (42) 

^„(:c) =AA„sech«-^/2(a;)e^''^'^^'^**"("'"'^(^')p^-''^-9'^^-«)(isinh(a;)). (43) 
Therefore, we define the complex Dirac-Scarf four vector potential as: 

{eV{r), eA{r)f) = ^0, (Ctanh(r) - {riR + zr//)sech(r) - -)f^ . (44) 

Here, (, rjR, rji e K. By making use of (18), we find the effective complex Dirac-Scarf potential. 



{VR + ^VIf-{C+^/2r + l/i) 



sech^ (r) - 2(C + 1 /2) (t/h + «?7/)sech(r) tanh(r) + ■ 

(45) 



7 



By substituting (45 ) in (18), we have: 
dr'^ 



'^^ + [(r?ij + %)'-(C + l/2)' + l/4)]sech2(r) 



- 2(C + l/2){r]R + i?7/)sech(r) tanh(r) 



4>l{r) = {Et^-m'-CW^{r). (46) 



Comparing (46) with Schrodingcr equation for complex Scarf potential, (40), gives the following 
real relativistic energy spectrum for complex Dirac-Scarf potential as: 

El^-m^-e^-{C-nf, (47) 

and the corresponding eigenfunction as: 

V)„(r) = A/;sech^ (r)e'''^ arctan(sinh(r))p^-,7/-9,'7/-9) (j sinh(r)). (48) 

Therefore, for > the real relativistic bound state energy eigenvalue is: 



^r„ = \/m2 + C2-(C-n)2), n = 0,l,2,...,n„„^, (49) 

where, Umax is the largest positive integer which satisfy 



CI < (50) 



6 The complex Dirac-Poschl- Teller Potential 

The complex Poschl- Teller Potential is [22] 

= M{M 1) _ NiN^l) 

smh^{t) cosh2(t) ^ ' 

where t = x - ie, x G (-oo,+oc), e e (0,7r/2) and M,N G 5R. In this case, V'-'^'^{x) is not 
singular at x = 0. In [22], the real eigenvalues and corresponding eigenfunctions of the potential 
(51) were found by solving the Schrodinger equation. It is found 

En = -{2n+aN+TM+{T-a)/2f, n = 1, 2, Umax < -i(aiV+rM+(r-a)/2), 

(52) 

where t = a = ±1 and 

^„(a;) = Ar„ sinh^^(t) cosh^^+i(i)P^^^-V2.<^JV+i/2) [cosh(2t)]. (53) 

So, by defining the complex four vector as: 

{eV{r) , eA{r)f) = (o, (C tanh(i) - j? coth(i) -'j)f^, (54) 

where t = r — ie,r G [0, oo), e G (0, 7r/2) and ( and i] are real. We can obtain Vfjj^'^ as: 

K/r^W = -C(C + l)sech\t) + 7,(7, + l)cosch2(i) + (C - rjf. (55) 
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Therefore, the second order differential equation for the upper spinor component of Dirac 
equation for Complex Dirac-Poschl- Teller is: 

(f ??(?? + 1) C(C + 1) 2 ^2l n (K(i\ 

"^ + ^1^-^^-^''™-"^ -(C-.)J0„W=O. (56) 

By comparing(56) with the associated Sclirociinger equation of the complex Poschl- Teller po- 
tential, (51), we can obtain the relativistic real energy spectrum as : 



Em = Vm^ - (C - V)^ - (2n + + ttj + (r - a)/2)^, (57) 
where n = 1, 2, Umax, which, rimax satisfy 



2nmax < ~ (C - - (f^C + r?7 + (r - a)/2), (58) 

and also the upper spinor component, (t)n{r) , is: 

,^;^(r) =Ar„sinh^''(i)cosh'^«+^(i)P^^''-^/2,aC+i/2)[pogl^(2t)], (59) 
where AT^ is normalization constant and Pn'^'"^ is Jacobi polynomial. 

7 Conclusion 

The non-Hcrniitian quantum theories have been studied for many complex potentials. It is 
observed that a relativistic quantum theory need not to be Hermitian to have a real spectrum. 
In this paper we obtain the real relativistic energy eigenvalues of a Dirac charged particle 
in complex static and spherically symmetric potentials. We show that these complex Dirac 
potentials have exact solution for all value of k{k is angular momentum quantum number). 
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